Geometrically induced singular behavior of entanglement 
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We show that the geometry of the set of quantum states plays a crucial role in the behavior of 
entanglement in different physical systems. More specifically it is shown that singular points at the 
border of the set of unentangied states appear as singularities in the dynamics of entanglement of 
smoothly varying quantum states. We illustrate this result by implementing a photonic parametric 
down conversion experiment. Moreover, this effect is connected to recently discovered singularities 
in condensed matter models. 
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Entanglement, a genuine quantum correlation, plays a 
crucial role in different physical situations ranging from 
information processing [jj to quantum many-particle 
phenomena [2J. Similarly to thermodynamics, smooth 
variations of controllable parameters which characterize 
a physical system may lead to singular behavior of entan- 
glement quantifiers. In some cases, in similarity to quan- 
tum phases transitions Q , these singularities are attested 
by abrupt changes in the quantum state describing the 
system. However, unexpected singularities may appear 
even when the quantum state varies smoothly [4(. Here 
we demonstrate how the geometry of the set of unentan- 
gied states can be related to singular behavior in physical 
phenomena. In particular we show that singularities at 
the boundary of this set can be detected by measuring 
the amount of entanglement of smoothly varying quan- 
tum states. 

Entangled states are defined as the states of composed 
quantum systems which cannot be written as a convex 
sum of products of density matrices for each composing 
part [f|. Separable states, on the other hand, admit such 
a representation and form a convex, closed set with posi- 
tive volume (for finite dimensional systems) [6|. This set, 
henceforth designated by S, is a subset of D, the set of 
all density matrices (S C D), which is also convex and 
closed. Therefore, a natural geometric way to quantify 
entanglement is to see how far - using some definition of 
distance on the state space - an entangled state is from 
the set S. This has been carried over for a variety of 
notions of distance, generating different measures of en- 
tanglement [7j . One of these geometric quantifiers is the 
random robustness, Rr, defined for any state p as the 
minimum s (s > 0) such that the state 
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motivation is clear: a represents a mixture of p with the 
random state 7r, and Rr (p) quantifies how much of this 
noise must be added to p in order to obtain a separable 
state. 

The main result of this Letter is to show that Rr can 
be used to investigate the shape of the boundary of S, 
dS. The principle is to take an entangled state depend- 
ing smoothly on one parameter q and compute Rr as a 
function of q. The one-parameter-dependent density ma- 
trices p(q) can be seen as a curve in the set of quantum 
states as shown in Fig. [TJ Singularities at dS will show 
up as singularities in Rr (p(q)). This statement is gen- 
eral for any finite dimension and will be formalized by 
the contrapositive: if dS is non-singular, then Rr (p(q)) 
is also non-singular. More precisely: 

Proposition 1. Let D be a closed, convex set. Let S C 
D also be closed and convex, with tt a point in the interior 
of S . IfdS is a C m manifold and the states p(q) describe 
a C m curve in D with no points in the interior of S 
and obeying the condition that the tangent vector p' (q) 
is never parallel to 7r — p(q), then RR(p(q)) is also a C m 
function. 

One must remember that a manifold is called C m if 
it can be parameterized by functions with continuous 
derivatives up to order m 



25] 



1 + s 



is separable (ir = I/d, where / is the identity matrix and 
d the total dimension of the state space) Q . The physical 



The reader can change 
C m by smooth, in the usual sense of C°° , with almost no 
loss (actually, we use smooth throughout this Letter in 
the less precise sense of "as regular as necessary"). Other 
topological remarks before the proof: the fact that S has 
interior points implies that S and D have the same di- 
mensionality (since there is an open ball of D contained 
in S), and the proof will use the notion of (topological) 
cone, which simply means the union of all segments from 
a given point V to each point of a given set A: this is 
called the cone of A with vertex V. 

Proof: The geometrical situation leads to the cone, given 
by (p,q) * P K + (1 — p) P (q), P S [0, 1]. The condition 
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on the tangent vector (together with the fact that it is 
interior to S, while p (q) has no point in this interior) is 
sufficient for this cone to be C m , except at the vertex w, 
at least locally in q. 

As S is bounded and convex, and it is in its interior, 
every straight line from it crosses dS exactly once. As 
p (q) has no point in the interior of S, this crossing always 
happens for < p < 1. Denote this crossing value by 
p c (q). The curve q i-> p c (q) n + (1 - p c (q)) p (q) is C m , 
implying p c is a C m function of q. 

The random robustness is given by Rji (p (q)) — ^ . 
As p c < 1, we also obtain that Rji is a C m function of" q. 
□ 

We insist on the interpretation: Proposition [1] means 
that any singularity in Rr for a well choosen path p (q) 
reflects singularities in dS. 




physical setup? We proceed to answer positively both 
questions by showing physical processes where a singu- 
larity in dS is revealed by monitoring the entanglement 
of a given system. 

First, let us consider a general system of four qubits a, 
b, A, and B, subject to the following Hamiltonian [l9j |: 

H = H aA + H bB , (3) 

where 

H» v = fa + fa + |(a^< + a$a1). (4) 

Here ct+ = (a x + ia y )/2 and er_ = (a x — ia y )/2, where 
o x ; cT y and a z are the usual Pauli matrices. This scenario 
can be realized in many systems, like cavity QED [20| . 
trapped ions [2l[, and quantum dots We set the 

initial state to be \tp(t = 0)) = \$+} ab <8> | V E'+) AB , where 
qubits ab are in the Bell state |$+) = (|00) + |ll))/v / 2 and 
qubits AB are in the orthogonal Bell state \^f+) = (|01) + 
|10))/\/2. Hamiltonian ([3]) induces a swapping process 
which leads (in the interaction picture) to the following 
temporal evolution for the subsystem AB, obtained by 
tracing out the subsystem ab: 



FIG. 1: State space. The dotted line represents the path 
p(q) followed by p when parameter q is changed. It is worth 
noting that S can present singular points in its shape and to 
remember that the "true" picture is much subtler, given the 
large dimensionality of even the simplest example 



PAB {t) = q (* + | + (l-q) |$+) ($+| , (5) 

where q — cos 2 (gt). For this state the negativity reads 
N{pA B (t)) - max{l - 2q,2q - 1} = |1 - 2q\ . (6) 



From this point on, we study the situation for two 
qubits, which is related to the performed experiment 
described here. In this case, Ref. [l(| shows that the 
Random Robustness is proportional to the Negativity 
{Af{p))[Xl], given by the absolute value of the sum of 
the negative eigenvalues of the partial transposed state. 
The negativity is a monotone under local operations and 
classical communication (l2l | and has the operational in- 
terpretation of a cost function under a certain class of 
operations 13L 14 1. 

At the same time, entanglement can be measured with 
the help of entanglement witnesses [l5[ ■ These are Her- 
mitian operators with positive mean value for all sep- 
arable states, but with a negative mean value for some 
entangled states [lf|. In fact, many geometrical entangle- 
ment quantifiers are directly related to witness operators 
17 1 . In the particular case of two qubits f26| . we have 
that for every entangled state p fic|. 



2N{p) = R R (p) = -2 mjn Tr(Wp), 
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where W is the set of entanglement witnesses W with 
TrW = 2. 

At this point we might ask some natural questions. 
Is there in fact any singularity in the shape of S? In 
the affirmative case, does this singularity appear in any 



This function presents a singularity for q = 0.5 (gt = 
n7r/4, with n odd) signaling then a singularity at dS. 

Another physical process which also produces the fam- 
ily of states ([5]) is the following simple quantum commu- 
nication task: Alice prepares a Bell state l^) and sends 
one qubit to Bob through a quantum channel; if this 
channel has a probability q of introducing a bit flip, and 
1 — q of no error at all, the state ([S]) is the output of the 



process [27 1 



To illustrate the dynamics given by Eq. ([S]), we have 
performed an optical experiment, shown in Fig. [2] In 
our experiment, twin photons maximally entangled in po- 
larization are generated in a non-linear crystal [23| and 
sent to an unbalanced Michelson interferometer, which 
is used to simulate the channel described above. The 
experiment works as follows: we produce a two-photon 
|* + ) state, send one of the photons directly to the de- 
tection stage, and the other to the (unbalanced) interfer- 
ometer. One of the arms of this interferometer does not 
change the polarization of the photon, and if the pho- 
ton went through this path the two photons would be 
detected in |^+). However if the photon went through 
the other path its polarization would be rotated in such a 
way that the final two-photon state would become |$+). 
We have made a tomographic characterization of the pho- 
tonic states corresponding to these two extremal points. 
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FIG. 2: Experimental setup: The state source is com- 
posed by a 2mm-thick BBO (/3-BaB2 04) nonlinear crystal 
(Ci) pumped by a cw krypton laser operating at 413nm, gen- 
erating photon pairs at 826nm by type II spontaneous para- 
metric down-conversion. Crystal Ci is cut and oriented to 
generate either one of the polarization entangled Bell states 
>]/_) or 1 4'+}. Walk-off and phase compensation is provided 
by the half-wave plate Ho followed by a lmm-thick BBO crys- 
tal (C2) together with two lmm-thick crystalline quartz 
plates (Z) inserted in one of the down-converted photon paths. 
The unconverted laser beam transmitted by crystal Ci is 
discarded by means of a dichroic mirror (U). The detection 
stages are composed by photon counting diode modules Di 
and D2, preceded by 8nm FWHM interference filters Fi and 
F2 centered at 825nm, and by circular apertures Ai of 1.6mm0 
and A2 of 3.0mm 0. Single and coincidence counts with 5ns 
resolving time are registered by a computer controlled elec- 
tronic module (CC). Polarization analyzers are composed by 
quarter- wave plates Qi and Q2, half- wave plates Hi and H2, 
followed by polarizing cubes Pi and P2. The State Source 
produces state |^-). For each pair, the photon emerging in 
the upper path goes straight to the polarization analyzer and 
to the detection stage 1. The lower path photon is directed 
by mirror M3 through the circular aperture A3 into the state 
mixer (an unbalanced Michelson interferometer), composed 
by the beam splitter BS, mirrors M4 and M5, quarter- wave 
plates Q4 and Q5, variable circular apertures A4 and A5, and 
by the half- wave plate H3, whose purpose is to compensate 
for an unwanted slight polarization rotation caused by the 
beam splitter. The quarter-wave plate Q4 is switched off 
which means that if the lower photon follows path labeled 
4, there is no change to its polarization and the half-wave 
plate H3 changes the state to 1^+}. On the other hand, if 
the lower photon follows path labeled 5, Q5 is oriented with 
the fast axis at 45° in order to flip its polarization. The path 
length difference, 130mm, is much larger than the coherence 
length of the down-converted fields, ensuring an incoherent 
recombination at BS. The pair detected by CC is in state 
+ (1 — q) I $+)(<!>+ 1 where q is defined by the rela- 
tive sizes of apertures A4 and A5. 



The reconstructed density matrices are displayed in Fig. 
[5J These two possibilities are then incoherently recom- 
bined, thus allowing the preparation of state Each 
preparation yields a different value for q with the corre- 




FIG. 3: (Color online) The reconstructed density matrices 
corresponding (ideally) to the states |$+) (A. real and B. 
imaginary parts) and 1^+) (C. real and D. imaginary parts). 
The attained fidelity for these states are, respectively, F$, = 
(<S>+\p |$+) « (92±3)%andF* + = (*+|p|*+) « (96±3)%. 



sponding optimal witness given by 



W opt = 



/-2|*+)($+| , for0<<7< 1/2, 
I- 2 ]*+)(*+ j , for l/2<q< 1. 



(7) 



For the family of generated states these two observables 
are the only candidates of optimal entanglement wit- 
nesses, so they are the only ones to be measured. In 
a more general situation, if less is known about the pre- 
pared state, much more candidate witnesses should be 
measured. The results are displayed in Fig. |U The blue 
curve in the figure shows the witnessed negativity mea- 
surement and its edge indicates the existence of singular- 
ities at dS. This experimental result shows the abrupt 
change in the optimal witness at the value q — i, which 
heralds the singularity in dS. As a proof of principles, 
each operator W is measured for the whole range of q, 
which yields the points bellow zero in Fig. 2] Note that 
the singularity occurs exactly for Rr = (q = 1/2). Ac- 
cording to our geometrical interpretation, this means the 
path followed by the parameterized state p(q) touches the 
border of S. This result must not be a surprise, since it is 
well known that in the tetrahedron generated by the Bell 
states (which we access in our experiment) the separable 
states form a inscribed octahedron (24|. 

The geometrical properties of entanglement discussed 
here give new insight into singularities found recently in 
the entanglement of condensed matter systems. Striking 
examples, dealing with entanglement properties of cer- 
tain spin-i models subjected to a transverse magnetic 
field h, are described in Refs. In these works, the 
two-qubit reduced state shows a singularity in entangle- 
ment for a particular field value h t far from the critical 
field of the respective model. As correlation functions, 
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FIG. 4: (Color online) Measurement of the mean value of 
both operators described in (0 for the full range < q < 1. 
Each W is expanded as a linear combination of products of 
local operators which are then measured independently. The 
blue continuous line corresponds to the theoretical value of 
H{p{q)) for the state p(q) = q |*+> <*+| + (1 - q) |$+> ($+[. 
Note that each W only witnesses entanglement for a restricted 
range of q values as predicted by the theory. T he local sin- 
gularity of dS is evidenced by the abrupt change of optimal 
W . Experimental errors are within the dots' sizes. 



ground state energy, and even reduced density matrices 
are all smooth at hf, there was no clear origin for these 
singularities. Our results offer an explanation by inter- 
preting the non-analyticities exhibited by entanglement 
as a consequence of geometric singularities at dS (28}. 

As previously mentioned, Rr can be used to probe dS 
in any finite dimensional system. For example, a previous 
work showed a singular behavior of Rr in three qubits 
systems Within the scope of our paper, we can 

interpret it as originated by a singularity at the border 
of the respective separable set. Note, however, that in 
this case, due to the higher dimensionality of the system, 
the singularity at dS occurs in the interior of D, with Rr 
showing a singularity at a positive value. 

To sum up, we have presented a method for probing the 
shape of the set of separable states. Singularities in this 
set were found and connected to non-analytical behavior 
of entanglement in different physical systems. It is an 
interesting open question to find physical implications of 
such singularities. 
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